Orbital order in Z11V2O4 
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In view of recent controversy regarding the orbital order in the frustrated spinel ZnV204 , we 
analyze the orbital and magnetic groundstate of this system within an ab initio density functional 
theory approach. While LDA+U calculations in the presence of a cooperative Jahn- Teller distor- 
tion stabilize an A- type staggered orbital order, the consideration of relativistic spin-orbit effects 
unquenches the orbital moment and leads to a uniform orbital order with a net magnetic moment 
close to the experimental one. Our results show that ab initio calculations are able to resolve the 
existing discrepancies in previous theories and that it is the spin-orbit coupling alongwith electronic 
correlations which play a significant role in determining the orbital structure in these materials. 

PACS numbers: 71.20.-b, 71.15.Mb, 71.70.-d, 71.70.Ej 



Correlated electronic systems involving transition 
metal oxides took the centrestage of condensed matter 
physics research for the last three decades because of 
their intriguing, often non-intuitive properties. Transi- 
tion metal spinel oxides with an additional complexity 
of a geometrically frustated lattice provide an exciting 
ground for the study of several competing interactions 
among spin, orbital and lattice degrees of freedom^ 2 . Be- 
sides being of fundamental interest, spinels have been also 
proposed for spintronics applications 3 -^. In the present 
work on ZnV204 we investigate the effect of competing 
spin, orbital and lattice degrees of freedom and show that 
density functional calculations provide an adequate and 
realistic ground to establish the dominant mechanism 
driving the orbital order in vanadium spinels. 

The orbital order in ZnV^CU as well as in other vana- 
dium spinels such as MgV204 and CdA^CU is presently a 
subject of considerable debate 5 . In order to understand 
the behavior of these compounds, various groups&££ 
have proposed alternative microscopic mechanisms which 
predict different orbital patterns. The ongoing debate 
has its origin in the complex nature of these systems with 
competing spin, orbital and lattice degrees of freedom. 
These systems have V 3+ ions in a spin 1 state charac- 
terized by double occupancy of the triply degenerate ti g 
(d xy , d XZl d yz ) orbitals. These partially filled t2 g orbitals 
leave the orbital degrees of freedom unfrozen opening up 
the possibility of orbital order. Moreover, the V-sites 
in the cubic spinel structure form a pyrochlore lattice, 
which gives rise to frustrated antiferromagnetic interac- 
tions among these sites 9 . In ZnV^CU the interplay of 
all these degrees of freedom leads to two successive phase 
transitions which involve structural, orbital and magnetic 
changes. At T$= 51 K, ZnV^C^ undergoes a structural 
phase transition where the symmetry is lowered from cu- 
bic to tetragonal with a compression of the V06 octahe- 
dron along the c axis^ and the system possibly orbital 
orders. The structural transition also lifts the geometri- 
cal frustration of the cubic phase making a way for the 
second transition at T/v= 40 K which is of magnetic na- 
ture and the system orders antiferromagneticall y 10 ! 11 . 



While the antiferromagnetic structure in ZnV^C^ at 
T < T^ has been unambiguously determined by neu- 
tron scattering experiments at low temperatures 12 , the 
nature of the orbital groundstate in the whole temper- 
ature range T < Ts is still unclear. Tsunetsugu and 
Motome 6 proposed that the Coulomb and exchange in- 
teraction between the magnetic ions as expressed in the 
Kugel-Khomskii Hamiltonian and the coupling to the 
Jahn- Teller lattice distortion are the driving mechanisms 
of the consecutive phase transitions. The groundstate or- 
bital order predicted by their model has an alternating 
singly occupied d xz and d yz orbital along the c direc- 
tion together with a singly occupied d xy orbital in all V 
sites (see Fig. [TJ right panel). They showed that this 
type of orbital order is consistent with the observed an- 
tiferromagnetic order at low temperature a 11 ! 12 as well as 
with the compressed V06 octahedron in the tetragonal 
phased. But this orbital order breaks the mirror re- 
flections in the planes (110) and (110) and the diamond 
glides in the planes (100) and (010) and hence is found 
to be inconsistent with the spatial symmetry IA\/amd 
predicted by X-ray scattering experiments on polycrys- 
talline sample a 1Q i 13 . Recent inelastic neutron scattering 
data on a single crystal sample 2 - seem to suggest though 
a breaking of the glide symmetry. These results have 
nevertheless not been further corroborated by other mea- 
surements. 

Tchernychyov 7 , in a separate theoretical work, consid- 
ering the fact that V 3+ has a non-negligible orbital mo- 
ment, proposed a single- ion model where the relativistic 
spin-orbit coupling plays the dominant role in driving 
the orbital order in the tetragonal phase. Including the 
Jahn- Teller effect as a perturbation in his model, Tch- 
ernyshyov predicted a uniform (ferro) orbital order with 
a complex orbital state d xz ±id yz at each V-site. This 
orbital order is consistent with the symmetry IA\/amd 
since it preserves both mirror and glide symmetries and 
hence is at odds with the prediction of Tsunetsugu et al 6 . 
Di Matteo et al. 8 in a subsequent work proposed a third 
alternative model where both the relativistic spin-orbit 
coupling and the Kugel-Khomskii model of spin-orbital 
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FIG. 1: (Color online) Schematic diagram showing the orbital 
order obtained within LDA+U in this work (left panel) and 
that proposed in Ref. (right panel). 
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FIG. 2: Total density of states for Zn, V, and O in both spin 
up and down channel in LSD A. 



superexchange are treated on an equal footing. They ob- 
served that for intermediate spin-orbit coupling strength 
the groundstate of the system shows a uniform orbital or- 
der compatible with Tchernychyov's prediction and with 
the symmetry IAi/amd. 

In view of the above mentioned disagreements among 
the various theoretical models, we analyze in this work 
the orbital and magnetic groundstate of ZnV2 04 in 
the tetragonal phase within an ab initio density func- 
tional theory (DFT) approach. Our calculations without 
spin-orbit coupling indicate the presence of a staggered 
antiferro-orbital order compatible with the spatial sym- 
metry of the tetragonal phase IAi/amd and we observe 
that correlation effects are essential for the description 
of this state. These results agree with the prediction of 
Tsunetsugu and Motome^ about the existence of a stag- 
gered orbital order except that we obtain a pattern which, 
contrary to Ref. 6, is compatible with the underlying 
I4i /amd symmetry (see Fig. [I]). Inclusion of relativistic 
spin-orbit effects in our DFT calculations changes though 
this picture considerably. For moderate to large values of 
the onsite electronic correlation U, a uniform orbital or- 
der analogous to the cubic phas o 14 i 15 -where the orbitals 
are tilted due to the trigonal distortion of the structure- 
is stabilized. Also, a large orbital moment antiparallel to 
the spin moment is obtained in the presence of strong on- 
site correlation which gives rise to a net magnetic moment 
in close agreement with the experimentally observed one. 
Finally, we confirm that the antiferromagnetic order ob- 
served in neutron scattering experiments 
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the ab initio DFT calculations. 

Our DFT calculations were performed within the local 
spin density approximation (LSDA) and the LSDA+U 16 
approach. The spin-orbit coupling was treated us- 
ing a scalar-relativistic basis and the second variational 
method 17 (LSDA+U+SO). All calculations were done 
with the full-potential linearized augmented plane-wave 
code WIEN2ki£. The atomic sphere radii were chosen 
to be 1.96, 1.99 and 1.78 a.u. for Zn, V and O respec- 
tively and 40 k points mesh in the irreducible wedge was 
considered for Brillouin zone integrations We used the 
structure data measured in Ref. |_10|. 



In Fig. [2]we present the spin polarized density of states 
(DOS) of Zn, V and O within the LSDA for both spin 
directions. The V t2 g up spin states are partially oc- 
cupied defining a metallic solution within the LSDA ap- 
proximation. The calculated magnetic moment per vana- 
dium is 1.48/iB- The above results show that V 3+ is in 
a high spin (5 = 1) state implying strong Hund's cou- 
pling (Jh) strength. Due to the c-axis compression of 
the V06 octahedra, there is a small but finite splitting 
of the d X y-\eve\ from the doubly degenerate d xz /d yz level 
(not shown here). 

ZnV204 is a Mott insulator, the inclusion of corre- 
lation effects are expected to improve on the incorrect 
metallic solution given by the LSDA approach. It is also 
well know n 19 i 2Q that in Mott insulators with orbital de- 
generacies as in our case, orbital order plays a crucial role 
in driving any collective Jahn- Teller distortion present in 
the system. Calculations within the LSDA+U approach 
-where orbital dependent potentials are included- should 
then be more appropriate to describe the gap and the 
possible orbital order in these systems. 

In Fig. [3] we present the projected DOS of the V 
^-states in the up spin channel calculated within the 
LSDA+U with U=b eV and J= 0.9 eV 21 . Note that the 
d xy up band now gets completely occupied with one elec- 
tron and the partially filled doubly degenerate d xz /d yz 
band observed in LSDA is split into two opening up a gap 
of 0.4 eV. The second electron now occupies the lower en- 
ergy up band which has both d xz and d yz character. In 
Fig. [6] (left panel) we show the corresponding bandstruc- 
ture. The orbital ground state described below is found 
to be stable with U values in the range 4-6 eV. 

We present in Fig. 0] a three-dimensional plot of the 
spin density of ZnV204 calculated within LSDA+U in 
the range of energies between -2 eV and the Fermi level. 
This range of energies includes V ti g and O-p states. A 
close inspection of Fig. [H reveals that all V ions in the 
tetragonal spinel structure have one filled t2 g orbital in 
the afr-plane ( d xy in the reference frame denned with 
the a, b axes directed along the basal V-O^ bonds and 
c axis directed along the apical V-O a bond). This fea- 
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FIG. 3: Projected V tig DOS in the spin up channel within 
the LSDA+U approach. 




FIG. 4: (Color online) Three dimensional electron density 
plots showing the staggered orbital order in this system. The 
inset shows the occupation of vanadium X>2 g orbit als in an 
elementary V4O4 cube. The isovalue used here is 0.6e/A 3 . 

ture is already observed in the projected DOS in Fig. 03 
The second vanadium electron occupies a t^g orbital with 
symmetry d xz + d yz and d xz — d yz in alternate afr-planes 
respectively along the c-direction (please note the four 
lobes in planes perpendicular to ab in Fig. 0]). The V4O4 
cube (Fig. [4] (inset)) shows two V with the second t2 g 
electron in a d xz + d yz symmetry and two V with the 
second ti g electron in a d xz — d yz symmetry. 

This orbital structure is therefore antiferro-orbitally 
ordered along c (A-type) and ferro-orbitally ordered in 
the ab plane. In comparison, in the high temperature 
cubic phase the t2 9 orbitals are split into singlet a\ Q and 
doublet e g due to a small trigonal distortion present in 
the system. The higher energy doublet is equally occu- 
pied at each V-site implying the absence of any orbital 
order unlike the tetragonal phase. Note that the trigonal 
distortion is also present in the low temperature tetrag- 
onal phase alongwith the tetragonal distortion. Here we 
discuss our results in the (d xy , d yzi d^) -basis instead of 
dig and e g ' . 

While the A-type orbital ordering observed by Motome 
and Tsunetsugu 6 for Z11V2O4 with alternately occupied 
d xz and d yz orbital along the c-direction (Fig. [I] right 
panel) is found to be inconsistent with the underlying 




FIG. 5: (Color online) Three dimensional electron density 
plots showing the uniform orbital order within LSDA+U+SO. 



space group symmetry of the crystal (I4i/amd), the or- 
bital order described in the present work is compatible 
with the crystal symmetry as it is obtained from a self- 
consistent field calculation on the crystal structure de- 
scribed by the IAi/amd spatial group (Fig. [I] left panel). 

The orbital order described above is driven by the com- 
bination of the staggered distortion present in the exper- 
imental lattice structure with correlation effects. First of 
all, the cubic to tetragonal transition with a compression 
of the V06 octahedron along the c-axis splits the triply 
degenerate \,2g states into a lower energy d xy -oih\tdl and 
higher energy doublet which are linear combinations of 
d xz and d yz orbitals (i.e. d xz ±d yz ), then the degeneracy 
in the higher energy doublet is lifted at each vanadium 
site due to a co-operative Jahn- Teller (JT)-like distortion 
(combination of trigonal distortion of the V sites) present 
in the experimental lattice structure which makes the O- 
V-0 angles to alternate between 85° and 95° in a stag- 
gered fashion along c-axis. Hence the faces of V4O4 cubes 
in the afr-plane alternately expand or compress along the 
c-axis (see Fig. 2]). LDA calculations, due to its known 
inabilities to describe this type of electronically driven 
cooperative JT distortio n 19 ^ fail to lift the degeneracy 
of the higher energy doublet and both d xz ± d yz orbitals 
are partially filled at each V-site. Consideration of U 
completely lifts this degeneracy and an orbital ordering 
is obtained as shown in Fig. [H 

In the next step of our calculations we perform a rel- 
ativists calculation including the spin-orbit interaction 
(LSDA+U+SO) in the second variational method. The 
magnetization direction is set along the c-axis as the 
experimentally observed spontaneous magnetic moment 
points to this direction and also this is the high symmetry 
direction of the tetragonal crystal structure. Within the 
LSDA+U+SO approach a uniform orbital order (see Fig. 
[5|) is favoured in contrast to the staggered A-type order 
found within the LSDA+U only approach. The orbital 
state at each site results from the mixture of the d xz -\-id yz 
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FIG. 6: Bandstructure in the energy interval (—3, 3)eV within 
LSDA+U (spin up channel) (left) and LSDA+U+SO (right). 
The bands are shown along the path Z(0, 0, 2tv/c)— G(0, 0, 0) — 
X(ir/a, 7r/a, 0) — N(7r/a, 0, tt/c) — Z. Note the narrowing of 
the broad V d unoccupied band ( energy between 0.4 eV and 
1.8 eV) when considering SO. 

and d xz — id yz orbitals that were found to be staggered 
at alternate sites within the LSDA+U. Due to the partial 
occupation of both d xz + id yz and d xz — id yz orbitals, the 
system has a finite \ z = 0.75 \±b moment. Since this or- 
bital magnetic moment is antiferromagnetically coupled 
to the spin only magnetic moment (1.69 hb) we obtain 
a total magnetic moment of 0.94/i^ per vanadium ion. 
We note that the good agreement with the experimen- 
tally observed magnetic moment (0.63 /is) 10 implies an 
important contribution of the orbital degrees of freedom 
in the present system. In fact, LSDA+SO only calcula- 
tion provides a small orbital moment of about 0.04 \ib 
due to the inability of LSDA to treat orbital polarization 
effects 22 . Including atomic orbital correlations on top of 
LSDA+SO scheme through the Hubbard U allows for an 
appropriate description of the atomic orbital state. 

In order to understand the effects of the spin-orbit cou- 
pling, we compare in Fig. [6] the bandstructure within 
LSDA+U (left panel) and LSDA+U+SO (right panel). 
In case of LSDA+U the bands are shown in the spin 
up channel because the bandgap is formed in this chan- 
nel whereas the spin down bands are about 1 eV above 



the Fermi level. The significant increase in the bandgap 
within LSDA+U+SO approach (1.0 eV ) can be un- 
derstood by the fact that the small but finite orbital 
moment (~ 0.04/i#) obtained within LSDA+SO is driven 
to about 0.75 hb by the orbital polarization effects of 
U. With this large orbital moment, the spin-orbit inter- 
action term in the Hamiltonian pushes the unoccupied 
d-bands above the Fermi level further up in energy. We 
also observe a decrease in the bandwidth of these bands 
compared to the case of LSDA+U and thereby a sub- 
stantial increase in the bandgap. 

An experimental measurement of the band gap in this 
system could further establish the role of spin-orbit in- 
teraction and atomic orbital correlation and thereby cor- 
roborate the type of orbital ordering. Our calculations 
show that this (uniform) orbital ordering persists in the 
antiferromagnetic phase observed in INS 1Q i n . 

In summary, we have investigated the orbital order in 
3d 2 magnetically frustrated systems by performing ab- 
initio DFT calculations on ZnV204 and observe that al- 
though strong onsite electronic correlation effects in the 
presence of co-operative JT distortion stabilize an A- 
type staggered orbital order, the consideration of rela- 
tivists spin-orbit effects unquenches the orbital moment 
and leads to a uniform orbital order state with a large 
orbital moment. This orbital moment couples antiferro- 
magnetically to the spin moment of V with a net mag- 
netic moment close to the experimentally observed one. 
This study shows that for a realistic description of the or- 
bital order in ZnV204 it is the interplay between electron 
correlation, spin-orbit coupling and co-operative JT dis- 
tortions what drives the orbital order. Our results help 
to resolve the differences among all the previous theo- 
ries and present a consensus picture of the orbital order 
in ZnV204 compatible with the experimental findings. 
Furthermore, this work shows direct evidence from ab- 
initio calculations of the importance of spin-orbit effects 
in 3d 2 systems. 
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